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ABSTRACT 

We predict that there is a population of low-luminosity dwarf galaxies orbiting within the halo of the Milky Way 
that have surface brightnesses low enough to have escaped detection in star-count surveys. The overall count of 
stealth galaxies is sensitive to the presence (or lack) of a low-mass threshold in galaxy formation. These systems have 
luminosities and stellar velocity dispersions that are similar to those of known ultrafaint dwarf galaxies but they have 
more extended stellar distributions (half light radii greater than about 100 pc) because they inhabit dark subhalos 
that are slightly less massive than their higher surface brightness counterparts. As a result, the typical peak surface 
brightness is fainter than 30 mag per square arcsec. One implication is that the inferred common mass scale for Milky 
Way dwarfs may be an artifact of selection bias. If there is no sharp threshold in galaxy formation at low halo mass, 
then ultrafaint galaxies like Segue 1 represent the high-mass, early forming tail of a much larger population of objects 
that could number in the hundreds and have typical peak circular velocities of about 8 kms~^ and masses within 300 
pc of about 5 million solar masses. Alternatively, if we impose a low-mass threshold in galaxy formation in order to 
explain the unexpectedly high densities of the ultrafaint dwarfs, then we expect only a handful of stealth galaxies in 
the halo of the Milky Way. A complete census of these objects will require deeper sky surveys, 30m-class follow-up 
telescopes, and more refined methods to identify extended, self-bound groupings of stars in the halo. 

Subject headings: cosmology: theory — dark matter — galaxies: formation — galaxies: halos — 
methods: iV-body simulations 



1. INTRODUCTION 



Approximately twenty-five new dwarf galaxy com- 
panions of the Milky Way (MW) and M31 have been 
discovered since 2004, more than doubling the known 
satellite population in the Local Group in five years 
(IWillman et aLl [20051: IZ^ker et al.l[2006l: lGrillmaidl2006l: 



Maiewski et al-l 120071: iBelokurov et all 120071 iGrillmairl 



20091: IBelokurovet al.l [20091: IMartin et al.l 120091) ? The 

majority of these newly-discovered dwarfs are less lu- 
minous than any galaxy previously known. The most 
extreme of these, the ultrafaint MW dwarfs, have lumi- 
nosities smaller than an average globular cluster Ly — 
10^ — lO"* Lq, and were discovered by searches for stellar 
overdensities in the wide-field maps of the Sloan Digital 
Sky Survey (SDSS) and the Sloan Extension for Galac- 
tic Understanding and Exploration (SEGUE). Follow- 
up kinematic observations showed that these tiny galax- 
ies have surprisingly high stellar velocity d ispersions for 
their luminosities and sizes ia^ ~ 5 km s~^. [ Martin et al.l 
[2007t ISimon fc Gehal[2007t IGeha et all[200l " and subse- 
quent mass modeling has shown that the y are the most 
dark matter dominate d galaxies known (IStrigari et al.l 
120081: iWolf et al.ll2010l ). Remarkably, these extreme sys- 
tems are not only the faintest, most dark matter domi- 
nated galaxies in the universe but they are also the most 
metal poor stella r systems yet studied ([Kirbv et al.ll2008l : 
IGeha et al.l[2009l) . 

Perhaps the most exciting aspect of these recent dis- 
coveries is that they point to a much larger popula- 
tion. Detection in the SDSS is co mplete only to ^ 50 
kpc for the l east l uminous dwarfs ([Koposov et al.|[2008[ : 
[Walsh et al.l [2009() and straightforward cosmologically 
motivated luminosity bias and coverage corrections sug- 

^ Center for Cosmology, Department of Physics and Astron- 
omy, The University of California at Irvine, Irvine, CA, 92697, 
USA; bullockOuci.edu 



gest that there are between ^ 200 and ~ 500 ultra- 
faint dwarf galaxies within 400 kpc of the Milky Way 
([ToUerud et al.[ [2001 . with an ever increasing number 
beyond. Importantly, the luminosity-distance detection 
limits only apply for systems with peak surface bright- 
ness obeying /^y < 30 mag arcsec"^ ([Koposov et al.l 
I2008D . Any satellite galaxy with a luminosity of L ~ 
1000 Lq and a projected half-light radius Ra larger than 
about 100 pc would have evaded detection with current 
star-count techniques regardless of its distance from the 
Sun. 

Here we argue that there is likely a population of dwarf 
galaxies surrounding the Milky Way (and by extension, 
throughout the universe) that are so diffuse in stellar den- 
sity that they would have thus far avoided discovery. Our 
predictions rely on the fact that the effective radius i?e 
of a dark matter dominated, dispersion-supported galaxy 
with fixed stellar velocity dispersion will increase as its 
dark matter halo mass decreases. One implication of 
this idea is that the known ultrafaint dwarf spheroidals 
(dSphs) may represent the high (dark matter) mass tail 
of a larger distribution of stealth galaxies. These undis- 
covered systems should preferentially inhabit the small- 
est dark matter subhalos that host stars (with maximum 
circular velocities l^nax 15 kms~^) and their pos- 
sible presence should be accounted for in models that 
attempt to understand the satellite-subhalo problem in 
Cold Dark Matter ( CDM) models fKlypin c raL|[T999l: 
Moore et al.' '1999^ 'Bullock et al.' 2000; Stochr et aU 
2002: Zcntncr & Bullock 2003; Kravtsov ct al. 2004^ 
Strigari et al. 2007b|; iMaccio et al.l 12 009: B usha et alj 
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ravtsov,20ialOkamoto fc Frenk,2009l ). 



The link between subhalo mass and galaxy surface 
brightness may be particularly important for the quest 
to identify 'fossils ' of the reionization epoch in the lo- 
cal universe (e.g. [Ricotti et al.[ [2002t [Ricotti fc GnedirJ 



20051 IGnedin fc Kravtsovl [20M lOrban et all [20081: 



Ricotti et all l2008t iMadau et al.lfeoOSHBovill fc Ricottil 



20091: iMufioz et al l 120091: iRicottil l2oTol ). The smallest 
halos {Vynax ^ 15 kms"^) are the ones that experi- 
ence suppression after reionization. Before reionization, 
stars can form in these systems via H2 cooling provided 
sufficient H2 is available. Conclusive identification of 
these H2-cooling fossils in the Local Group would pro- 
vide present-day laboratories for studying first-light star 
formation. Unfortunately, according to our estimates be- 
low, these low-mass, first-light fossils will likely be too 
diffuse to discov e r read ily in star-count surveys. Indeed, 
iBovill &: Ricottil (|2009| ) have shown that low-luminosity 
fossil dSph galaxies produced in self-consistent cosmolog- 
ical simulations do tend to have very low surface bright- 
ness, such that they would not have been discovered in 
current surveys. 

Stellar kinematic samples in dSph satellte galaxies 
provide a means to directly constrain total dark mat- 
ter halo masses or Knax values on an object by ob- 
ject basis. Unfortunately, this task is not straight- 
forward because dSph dark matter halos are expected 
to extend well beyond their stellar radii. Neverthe- 
less, by constraining the mass of a galaxy within its 
stellar extent, one can estimate (for example) a Vmax 
value by imposing a prior assumption that their density 
profiles b ehave as p redicted for subhalos in C DM sim- 
ulations (IZentner fc B ul lock 120031: iHavashi et al 200l 
iKazantzidis et all 12004 iPefiarrubia et al.l l2008al ). Re- 
markably, these studies indicate that all of the known 
MW dSph satellite galaxies with well-studied kinematic 
samples are embedded within massive dark matter sub- 
halos (Knax ^, 15 kms~ ^: Wolf et al., in preparation; 
IPefiarrubia et al.l l2008aD . and this includes ultrafaint 
dwarfs. Keep in mind, however, that these determi- 
nations rely on significant extrapolations: the rotation 
curve of a Knax = 15 kms~^ subhalo typically peaks at a 
radius of ~ 1700 pc (Springel et al. 2008; Diemand et al. 
2008). The median three-dimensional half-light radius 
for the sample of dSphs we consider here (see Wolf et al. 
2010, Table 1) is about 300 pc, therefore the extrapola- 
tion is fairly large. 

Another way to study dSph halo potential well depths 
is to consider masses within central radii that are simi- 
lar in extent to the stars in the galaxies (jStrigari et al.l 
I2007a^ but still large enough that they can be re- 
solved directly in numerical simulations (jStrigari et al.l 
I2007U 120081 ). Today's state-of-the art N-body simula- 
tions cannot resolve subhalo densities below about 300 
pc (Springel et al. 2008; Diemand et al. 2008). As a for- 
tunate coincidence, 300 pc is also the median half-light 
radius for the popu lation of MW dSphs with well-studie d 
kinematic samples (jStrigari et al.l2008HWolf et al.ll2010l). 
For th ese practical (not physical) reasons Istrigari et al.l 
(j2008l ) studied the integrated mass within 300 pc (M300) 
for all of the MW dSphs with large stellar kinematic sam- 
ples. They found A/300 — lO^A/© and no evidence for 
a relationship between halo M^qq and total luminosity: 
M300 oc 2,o.03±o.03_ j'Yie implication is that the tiniest 
satellite galaxies have dark matter densities indicative of 
dark matter halos that are at least as massive as those 
of systems 10,000 times more luminous. 

The normalization of the Strigari et al. relation at 



M300 — IO^Mq is indicative of the central densities of 
massive dark matter subhalos (I^nax ^ 15 km s^^). As al- 
luded to above, this ma ss-scale is fa irly e asy to explain in 
ACDM mod els (Strigari et alll20 08: Mac cio et aLl|2009[ 
Busha et al. 2009; Munoz et al. 2009; Kravtsov 20ld 
Qkamoto fc Frenk ,2009.: Stringer et aD 
the lack of observed correlation be- 



Li et al.. .2009.: 
20091). However, 



tween L and mass is quite unexpected. To put the lack 
of measured slope in the M300 — L relation in perspec- 
tive, consider the relationship between dark matter halo 
Vmax and galaxy luminosity required to match the faint- 
end slope of th e galaxy lurninosity function: L cx V^,^iax 
with 6 = 7.1 (jBusha et al.ll2009f) . TuUy-Fisher stud- 
ies suggest a relation that is even l ess steep for more 
massive systems, with 6 4 (e.g., iStark e t al.l 120091 : 
ICourteau et 10120071 : iMcGaugh fc Wo lf 2010). For dark 
matter halos of interest, with maximum circular veloci- 
ties Knax — 15 — 45kms~^, w e expect M300 oc Vma x (as- 
suming NFW fits to halos in iSprmgel et al.l 120081) such 
that the observed trend M300 oc L"""^ would naively im- 
ply L oc V^ax with 6 ~ 33. This is a much steeper re- 
lationship than we expect from more luminous systems. 
The lack of inferred relationship between dSph luminos- 
ity and total halo mass is not an artifact of the specific 
choice of 300 pc for the mass comparison. If one performs 
the same comparison betwe en galaxies using a smaller 
benchmark radius (100 pc, iStrigari et al.l I2008D or us- 
ing t he 3d half-light radius for each galaxv ifWolf et ahl 
I2OIOI) or using the 3d mass within the 2d half-light ra- 
dius of each galaxy (j Walker et al.ll2009l ) then one reaches 
the same conclusion: there is no observed trend between 
inferred total halo mass (or Vmax) and luminosity. 

One possible explanation for this lack of mass trend 
is that it refiects a scale in galaxy formation, where 
the scatter in L at fixed Vnax becomes very large, as 
might possibly be explained by feedback due to photoion- 
ization or suppressi o n b elow the atomic coolin g limit 
(iRicotti et al.l I2OO8I: IStr igari et al. 2008; Macci 6^et al.l 
I2009t lOkamoto fc Frenfl 120091 : iStringer et al.ll200l and 
§3 below). Another possibility, outlined below, is that 
the lack of an observed trend between mass and lumi- 
nosity is the product of selection bias: most ultrafaint 
galaxies do inhabit halos with Af3oo < lO^M©, but they 
are too diffuse to have been discovered. 

In the next section, we explain why we expect surface 
brightness selection bias to limit the discovery of satel- 
lite galaxies in small subhalos. In §3 we use a simple 
model to estimate the number of low surface brightness 
stealth galaxies within the vicinity of the Milky Way. 
Our estimates rely on the public subhalo catalogs pro- 
vided by the Via L actca 2 (VL2) N-body simulation 
group (jPiemand et a l. .200^. We explore two models. 
Our Fiducial Scenario connects each subhalo's mass and 
accretion time to a galaxy luminosity L by extrapolat- 
ing the halo mass-light relationship required to match 
the asymptotic slop e of the galaxy stellar mass function 
(iMoster et al.ll2009j ). Our secondary model (Threshold 
Scenario) explores a scenario where galaxy formation is 

^ Note that there is an observed trend between the mass within 
the half-light radius of each galaxy and the galaxy's luminosity, 
but this simply reflects the fact that brighter galaxies have larger 
half-light radii. This trend is perfectly consistent with each of the 
galaxies be ing embedded with i n halos of approx imately the same 
total mass HWalker et al.|[2009l : IWolf et al.|[20Toh . 
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Fig. 1. — Projected helio-ccntric distance vs. V=band luminosity (left) and half-light radius Re vs. V-band luminosity (right) and for 
Milky Way dSph galaxies. The lines in the left panel show SDSS luminosity completeness limits from Walsh et al. (2009) and an estimate 
for LSST completeness foUowring ToUerud et al. (2008). In the right panel, galaxies above the solid line, with surface brightness fainter than 
/I = 30 mag arcsec^^ are currently undetectable. For reference, the dashed line in the right panel corresponds to = 35 mag arcsec"^. 



truncated sharply below a characteristic dark halo mass 
scale. We present our findings in §4 and conclude in §5. 

One of the goals of our model is to investigate how 
surface brightness detection limits can affect our inter- 
pretation of dwarf subhalo mass-luminosity trends. We 
will use M300 as our primary mass variable for charac- 
terizing subhalo masses. As discussed earlier, M^qq can 
be measured directly for the N-body subhalos, but re- 
quires less extrapolation than Knax when comparing to 
observational data. For the observational comparisons, 
we will use M300 values derived by Strigari et al. (2008) 
for MW dwarfs. We emphasize that the choice of 300 
pc is a practical one, with no special physical meaning, 
other than 300 pc is the median half-light radius of the 
dSph galaxies in our sample. Integrated masses within 
smaller radii cannot be resolved in the simulation. Eight 
of the nineteen galaxies we consider have 3d half-light 
radii smaller than 290 pc and there are four galaxies that 
do not have at least one kinematic stellar tracer beyond 
this radius. In these cases we are relying on the CDM- 
motivated prior of Strigari et al. (2008) to extrapolate 
masses out to 300 pc. Such an extrapolation is perfectly 
reasonable (and inevitable because of resolution) as long 
as our aim is to compare to predictions from the same 
theory that motivates the prior (as is the case here) . The 
assigned error bars on the measured M300 take into ac- 
count uncertainties in the extrapolation encompassed by 
the theory, including an allowance for exponential mass 
truncation because of tides. We note that the instanta- 
neous tidal radius of the closest dSph, Segue 1, is much 
larger than 300 pc (Geha et al. 2009). Moreover, by 
examining the orbits of subhalos in VL2, Rocha et al. 
(2010, in preparation) find that subhalos chosen to have 
radii and masses consistent with Segue 1, have past or- 
bital trajectories and measured tidal radii that are larger 



than 300 pc in the vast majority of cases. 



2. MOTIVATIONS 



The lef t panel of F i gure 1 shows the MW dSphs as tab- 
ulated in lWolf et al.l ()2010[ ) plotted in the plane of helio- 
centric distance vs. V-band luminosity. The solid line 
(labeled SDSS) illustrates the distance to which dwarfs of 
a given lum inosity can be detec ted in SDSS with 90% effi- 
ciency fro mlWalsh et al.l (120091) Similar resuUs were pre- 
sented by iKoEQSOverainiool) . The upper line shows 
the same limit adjusted up by scaling to th e limiting mag- 
nitude of the full co-ad ded LSST survey (jTollerud et al.l 
I2008tllvezic et al.ll2008[ ). Clearly, the known dwarf galax- 
ies cluster at the current completeness edge of the dia- 
gram, indicating a high likelihood fo r future discoveries 
jKo posov et al.l l2008tlTollerud et al.l[200l iWalsh et al.l 
I2009h . 

The distance-luminosity completeness l i mits p resented 
by Walsh et al. (2009) and lKoposov et al.l (|2008l) are only 
applicable for systems with s urface brightness bri ghter 
than fly = 30 mag arcsec^ -2 (jKoposov et al.ll2008L and 
G. Gilmore, M. Geha, and B. Willman, private commu- 
nications). Systems more diffuse than this limit cannot 
be detected in SDSS with current methods, no matter 
their helio-centric distance. This phenomenon is illus- 
trated qualitatively in the right panel of Figure 1 , which 
presents the same set of MW dSphs in the plane of Rc 
vs. L. The solid line shows a constant peak (central) 
surface brightness for a Plummer profile 
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0.036 L0pc~ 



(1) 



and corresponds to fiy — 30 mag arcsec"^ for solar 
absolute magnitude Mqv = 4.83. As in the distance- 
luminosity figure, the tendency for many of the fainter 
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Fig. 2. — Mass within 300 pc vs. luminosity (left) and measured stellar velocity velocity dispersion, ct*, vs. luminosity (right) for Milky 
Way dwarf spheroidal galaxies (data points with error bars). In the left panel, galaxies within in the shaded region (below dashed line) will 
remain hidden {fj, > 30 mag arcsec"^) if they have > 4kms~^ ( o"* > 5kms~^). In the right panel, galaxies in the shaded region will be 
stealth (/i > 30) if they have Afsoo < 2 X IO^Mq and galaxies below the da shed hne will be s tealth if they have M300 < IO^Mq. The masses 
and velocity dispersions quoted here are taken from lStrieari et al.l II2008I ) and IWolf et al.l II2010I ). Since the time of those publications, it 
has come to our attention that the velocity dispersion (and hence mass) errors on Hercules (at L ~ 3 X lO^LfT^ . Eden et al.ll2009l ^ and Wl 
(at L ~ 10"^ -Lq, B. Willman & M. Geha, private communications) are likely underestimated because of membership uncertainties. 



dwarfs to "pile up" near the surface brightness detec- 
tion hmit is suggestive. There is nothing ruhng out the 
presence of a larger population of more extended sys- 
tems that remain undetected because of their low surface 
brightness. 

If a large number of diffuse, undetected galaxies do 
exist, they are likely associated with low-mass dark mat- 
ter subhalos. One can understand this expectation by 
considering an spherically-symmetric galaxy that is in 
equilibrium with stellar density distribution p*(r) and 
radial velocity profile ar{r) that is embedded within a 
gravitationally-dominant dark matter halo mass profile 
M{r). The Jeans equation is conveniently written as 



MM = ^(7. + 7. -2/3), 



(2) 



where /3{r) = 1 — cr^/cr^ characterizes the tangential 
velocity dispersion and 7* = — dlnp^^/dlnr and "fa- = 
— dlntJ^/dlnr. If we make the simplifying assumption 
that /3 = and ar(r) ~ cr^ = constant, with 70- ^ 1 then 
M{r) = r al ji,. For a fixed velocity dispersion, a 
more spatially extended profile (smaller 7^,) requires a 
lower mass at fixed radius. 

The same basic expectation follows in a more general 
context from the recent work of Wolf et al. (2010), 
who showed that the total mass of a quasi-spherical 
dSph galaxy within its 3d half-light radius ri/2 — 
1.3 Ro may be determined accurately from the luminosity 
weighted line-of-sight velocity dispersion cr^, for general 



1/2) 



3G 



Mass determinations at 



/3: M(7 

larger and smaller radii require an extrapolation of the 
mass profile from that point, but given a theoretical pre- 
diction for the mass profile shape M(t) one can perform 
this extrapolation by simply normalizing at r = ri/2- 
It is useful to rewrite the Wolf et al. (2010) mass esti- 

^ Under the assumption that the observed stellar velocity dis- 
persion remains fairly fiat with projected radius, as is the case with 
all of the well-studied systems. 



mator in terms of the implied circular velocity at r-ij^'. 



1/2) 



\/3(T*. 



(3) 



Consider then a galaxy with velocity dispersion Ui, 
and luminosity L embedded within a gravitationally- 
dominant dark matter halo described by a circular veloc- 
ity curve that increases with radius as an approximate 
power law: Vc{r) = V300 (''/SOO pc)". Equation [3] im- 
plies ri/2 = 300 pc (v^cr^/Vsoo)^^"- As an example, for 
an NFW halo (jNavarro et al.lll997| ) with r,, > 300 pc 
we have a — 1/2 and ri/2 V^m ^ -^^300- Clearly, the 
galaxy becomes puffier as we decrease M300 or V300 ■ One 
implication is that if a galaxy has a stellar density that 
is just large enough to be detected, another galaxy with 
identical L and cr^, will be undetectable if it happens to 
reside within a slightly less massive halo. 

Figure 2 provides a more detailed exploration of the re- 
lationship between halo mass parameters (M300 or Knax) 
and associated dSph observables cr^, L, and Re- Points 
in the left panel of Figure 2 presen t M300 vs. L f o r MW 
dSph galaxies, with masses fr om iStrigari et ahl (|2008f ) 
and luminosities updated as in lWolf et al.l ( 2blCl[ ). The 
right panel sho ws a^, ys. L for th e same galaxies culled 
from Table 1 of iWolf et all (poTol ). 

The shaded bands and dashed lines in each panel of 
Figure 2 illustrate the way in which surface brightness in- 
completeness may affect these diagrams. In determining 
these regions we have assumed each dSph is dark-matter 
dominated, such that its gravitating mass profile pro- 
duces an NFW circular velocity curve Vc{r) = V^-p-w[f)- 
Given the NFW shape, the rotation curve is fully spec- 
ified by its peak value Knax and the radius where the 
peak occurs r^ax (e.g., iBullock et al.l |2001[ ). We as- 
sume for simplicity that subhalos of a given Vmax map 
in a one-to-one way to a rotation curve shape using 
fmax = 650 pc (Vmax/lOkms"^)-^ '^^, which is indicative of 
median subhalos in high-resolution N-body simulations 
(intermediate between the normalizations of Springel et 
al. 2008 and Diemand et al. 2008). With this as- 
sumption in place, given a halo mass variable (e.g.. 
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-^300 or T^nax), wc may determine the implied half- 
light radius Re ~ 0.75ri/2 associated with any a^, using 



10' 



NFw(^l/2 



VScT^, (Equation [3]). 



In the right panel of Figure 2, galaxies residing in the 
shaded region are unobservable if they sit within dark 
matter halos less massive than M300 = 2 x 10'' or (equiv- 
alently for our assumptions) with peak circular velocity 
smaller than T^nax — 24kms~^. Similarly, galaxies resid- 
ing above the dashed line are too diffuse to be detected if 
they have M300 < lO^A/© or Knax S 14kms~^. Galaxies 
need to have deep potential wells if they arc to remain ob- 
servable at low luminosity for a^, ~ 5kms~^. If there are 
low-luminosity galaxies with A/300 values smaller than 
~ IO^Mq they would remain hidden as long as they have 
stellar velocity dispersions comparable to those of the 
known dwarfs. 

A related set of limits in the M300 ~ L plane is de- 
picted in the left panel of Figure 2. Galaxies sitting in 
the shaded band will have i?o too large to be observ- 
able if they have cr* > 4kms~^. Slightly hotter galaxies, 
with a-i, > 5 kms~^ will be unobservable if they sit below 
the dashed line. As expected, the hotter the galaxy, the 
deeper the potential well needs to be in order keep the 
stars confined to an observable surface brightness. We 
see that galaxies with L ^ 10'^ L© residing in a halos 
less massive than M300 ~ 8 x 10^ M© will be too diffuse 
to be seen if they have cr* = 5kms~^. Note that for 
L > lO^L© the constraint on allowed M^qq values is flat 
with L because halos smaller than this value are kinemat- 
ically forbidden via Equation [3l Specifically, kinematic 
mass determinations demand Vmax ^ Vio'*- 

Implicit in the above discussion is the idea that a 
galaxy's cr* can be considered independently of its halo 
mass. Dynamically, the only constraint is that a^, < 
Vma-x^/V^ (Equation [S]). One is more inclined to sus- 
pect that (T^ in an ultrafaint dSph is governed by star 
formation and galaxy formation processes, with an ab- 
solute minimum set by the effective temperature of the 
star forming ISM. Even for a very cold primordial ef- 
fective ISM temperature, Tjsm ^ 300K, we expect 
o'isM ^ 2kms~^, and this ignores turbulent and mag- 
netic pressure terms. The vast majority of globular 
clusters have stellar ve locity dispersions larger than this 
(jPrvor fc Mevlanlll993l ). Moreover, dark matter halos of 
all masses are expected to have ex perienced significan t 
mergers in their early histories (e.g.. lStewart et aLll2008D . 
These mergers would have heated (the oldest) stars be- 
yond any primordial pressurized motions, and this effect 
is indeed se en in cosmolo gical simulations of dwarf galaxy 
formation (|Ricottill2010l ). 

In the next section we consider the implications of a 
model where a^, is correlated with luminosity L in a way 
that tracks the observed relationship (right panel of Fig- 
ure 2). In principle, there could be a floor in the cr^, values 
allowed for dwarf galaxies. We do not impose such a floor 
in our calculations, but if one does exist, then our esti- 

* In deriving these regions, wc have explicitly assumed that the 
stellar systems are dark-matter dominated within their half-light 
radii. The same arguments cannot be applied to globular cluster 
systems, some of which do inhabit the shaded regions in the right 
panel of Figure 2 without any discernible dark matter halo. These 
systems have large velocity dispersions simply because they have 
very high stellar densities. 
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Fig . 3. — Model stellar mass - halo mass relation IjMoster et al.l 
12009) shown at three example redshifts. Our Fiducial Scenario 
assumes that the Mhaio"-^^* relation extrapolates smoothly to very 
small masses. Our Threshold Scenario imposes a sharp truncation 
mass at Mhaio = 5 X 10* (vertical dotted line) below which all 
halos are assumed to form no stars. 

mates could under-predict the fraction of stealth galaxies 
at low luminosities. 

3. MODEL 

We rely on the publicly released subhalo catalogs of the 
Via Lacte a II N-body simulatio n (VL2 hereafter) as de- 
scribed in lDiemand et al.l ()2008D . The sim ulation adopts 
cosm ological parameters from WMAP3 (jSpergel et al.l 
l2007f l and tracks the formation of a Milky Way size dark 
matter halo with a highest particle-mass resolution of 
4, lOOM© and force resolution of 40 pc. The main halo 
has a radius of 402 kpc, defined to enclose a mean density 
that is 200 times the mean matter density of the universe, 
and an associated mass of A/haio = 1-9 x lO^^Af©. The 
public subhalo catalogs include M^ooi Knax, and rmax pa- 
rameters for each bound system, as well as merger his- 
tory information that allows us to track the redshift of 
infall Zinf for each subhalo and to determine its maxi- 
mum attained mass Mmax prior to infall. Tests by the 
VL2 team suggest that the measured A/300 masses are 
good to about 20% (random) owing to resolution effects 
(J. Diemand, private communication). 

We assign light to each of our accreted dark matter 
subhalos by assuming that at each redshift z there 
is a monotonic relationship between halo mass A/haio 
and galaxy stellar mass A/*. This general approach 
is motivated by its demonstrated success in produc- 
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Of course, cosmological abundance matching cannot be 
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applied directly at the smallest stellar masses because 
of completeness issues. In our Fiducial Scenario we 
simply adopt the asymptotic A/haio - relationship 
suggested by the most complete stellar mass functions, 
which effectively assumes that there is no new (or 
abrupt) mass scale that truncates galaxy formation in 
small halos. We also explore a Threshold Scenario that 
imposes such a truncation scale (see below). 

For our Fiducial Scenario, we assign Mi, to each sub- 
halo by ex trapolating the fitting formula presented by 
iMoster et~a l. (2009) to small stellar masses. Moster et al. 
(2009) derived the relationship using N-body halo cata- 
logs together with observationally inferred stellar mass 
functions for M^, > lO^~^Af0 galaxies out to redshift 
z ~ 3. The implied (extrapolated) relationship between 
stellar mass and dark halo mass is presented in Figure 
3 for three example redshifts. We see that Mi, must de- 
crease at high redshift for a fixed Mhaio in order to ex- 
plain the evolving stellar mass function. Low-mass halos 
at high redshift have not had time to form as many stars 
as their z ~ counterparts. For our Threshold Scenario 
we adopt the same mapping for massive halos but we 
impose a sharp truncation in the Mi, - Mhaio relation at 
Mhaio = 5 X IO^Mq (dotted line in Figure 3). 

We assume that star formation is quenched in each 
subhalo at a time Tq after the redshift of accretion into 
the VL2 host. Specifically, subhalo light content is deter- 
mined at redshift Zq (set at a time Tq after the accretion 
redshift) using the appropriate Moster et al. (2009) map- 
ping with M = Mmax, the maximum mass each subhalo 
progenitor obtained prior to infall. If the subhalo is ac- 
creted at a time less than Tq before z = we adopt the 
Moster et al. (2009) relation at Zq = 0. For the fig- 
ures we present below we use a quenching timescale that 
is roughly a dynamical time for the host halo Tq = 2 
Gyr. We find that the value of Tq only affects our pre- 
dictions for the largest satellites Mi, > IQ^ Mq. For ex- 
ample, if we set Tq = 0, we under-predict the number 
of luminous satellites by a few, but the low-luminosity 
satellite count is largely unaffected. The main conclu- 
sions of this paper regarding the least luminous, stealth 
satellites are not sensitive to the choice of Tq. For all 
satellites, we convert from stellar mass to V-band lumi- 
nosity using AIi,/L = 2 [M(7) / L(7)] , which is typic al for 
Milky Way dSphs according to lMartin et all ()2008[ ) for a 
Kroupa IMF. 

In assigning luminosities to dark matter subhalos at 
the time of infall, we are effectively assuming that the 
majority of surviving dark matter halos in our model 
have lost very little stellar mass after infall (though al- 
most all of them lose dark matter mass). We do make 
a crude self-consistency check for stellar tidal mass loss 
below, but find that it is not significant in most cases. 
The idea that stellar mass loss has been minimal for 
dense, bound dark substructures in the halo is diffi- 
cult to test empirically with existing dwarf data (e.g., 
[Pcnarrubia et al. 2009). Therefore, in order to check 
whether this assumption is reasonable (at least in the 
context of the ACDM-b ased model we explore here ), we 
examined the output the lBullock fc JohnstonI (|2005l ) sim- 
ulations, which modeled hundreds of satellite accretions 
each for 11 ACDM merger histories using cosmologically- 
derived accretion times and orbits, set within an evolving 



Milky- Way-like disk and dark halo potential. We found 
that the overwhelming majority of satellite-subhalos that 
survived with bound dark matter halo cores (Vmax > 8 
km s~^) experienced no stellar mass loss. Indeed, 93% of 
the su rviving dwarf galaxies in the iBullock fc JohnstonI 
()2005D simulations lost less than 10% of their initial stel- 
lar material. The only surviving systems that show any 
stellar mass loss are those that have lost more than 90% 
of their initial dark matter mass. Stellar mass loss be- 
comes typical (affecting more than 50% of systems) only 
in the minority of subhalos that have lost more than 
95% of their initial dark matter mass. Moreover, we 
find no significant trend between current Galacto-centric 
distance and stellar mass loss in the simulated satellites. 
These facts provide encouraging support to the simple 
assumptions we adopt here. 

Once L is determined for each subhalo, we assign a 
stellar velocity dispersion by adopting the empirical re- 
lation shown in the right panel of Figure 2: 



ffi, ^ 6.9 kms 



L 



105 Lr. 



0.09 



(4) 



with a log normal scatter of A logj^g cr = 0.1 at fixed L (as 
measured in the data). Here, we are assigning stellar 
velocity dispersions to our subhalos at z = using a 
relation measured for real Milky Way dSphs. This means 
that the assignment is reasonable for subhalos, modulo 
the concern discussed above regarding stellar mass loss.0 

Once CT* is assigned, we determine ~ 0.75 r 1/2 us- 
ing VNFw{r 1/2) ^ VScFi, (Equation 2). The Kiax and 
''max values that define Vjv fw are those measured for 
each subhalo in the simulation (at the present day, not 
at the time of infall). Specifically, we are extrapolating 
the density profiles of the subhalos to radii < 300 pc that 
are not well resolved in the simulation. While this ex- 
trapolation is reasonable, it means that our derived Re 
values arc reliant on this assumption. 

For simplicity, wc assume that each dwarf galaxy fol- 
lows a Plummcr profile, with a peak surface density 
given by Equation [T] As discussed above, galaxies with 
Speak < 0.036 L0pc~^ are assumed to be undetectable 
with standard techniques. We note that if we impose a 
floor in allowed velocity dispersions near ct^, = 4 kms~^ 
then our results do not change dramatically. 

The final step in our procedure is a self-consistency 
check to see if our implied ri/2 values are small enough 
for the galaxies to be relatively unaffected by tidal strip- 
ping. In order to do this we estimate a tidal radius for 
each galaxy 7't and remove galaxies from our catalogs if 
fi/2 > ft on the assumption that most of their stars will 
have been tidally liberated (even though a dark mat- 
ter core remains bound). It is well known that a sub- 
halo's rotation curve should decline more rapidly than 
an NFW profile for r > r^ax because of tidal effects 

^ We do not account for any systematic surface brightness bias 
that would lead to high cr* systems being missed (as these are 
the systems that will have large i?c). By ignoring this effect we 
are systematically under- estimating the possible number of stealth 
galaxies. 

® Of course, if there is stellar mass loss, the least-bound (hottest) 
stars are stripped first, and this causes the veloci ty dispersion and 
luminosity to evolve in concert, as detailed in IPeiiarrubia et all 
l|2008bl ) for several example orbits. 
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Fig. 4. — Fiducial Scenario galaxy size vs. luminosity relation (colored points, loft); distance vs. luminosity relation (middle); and 
luminosity functions for different completeness cuts (right). Left Panel: The solid line corresponds to the current detection limit at a peak 
surface brightness of /.i = 30 mag arcsec"^ and the dashed line corresponds to = 35 mag arcsec"^, for reference. The small colored points 
represent model galaxies and the point-type scheme maps to the redshift of infall into the host dark halo: open triangles have ^jnf > 1.5, 
green diamonds have 0.5 < ^inf < 1.5, and blue plusses were accreted since z = 0.5. The black squares with error bars are known MW 
dSph galaxies. Middle Panel: The solid and dashed lines indicate luminosity-distance completeness in the SDSS and LSST, respectively, 
for systems with ^ < 30 mag arcsec"-^. The point types are the same as in the left panel. Right Panel: The symbols shown as X's reflect 
the current census of MW dSphs, corrected for the sky coverage completeness of SDSS as in ToUerud et al. (2008). The uncertainty in this 
correction corresponds roughly to the size of the symbols we use. The dotted line shows the predicted cumulative luminosity function of 
satellite galaxies that are bright enough to have been detected by SDSS according to the Walsh et al. (2009) limits. The dashed line shows 
the predicted luminosity function of all satellites with surface brightness meeting the < 30 mag arcsec"'^ threshold, most of which should 
be detectable by LSST. The solid line shows all satellite galaxies, including the stealth population. We see that the majority of ultrafaint 
dwarfs are expected to have surface brightnesses so low that they will avoid detection without revised techniques for discovery. 



(jKazantzidis et all |2004[) . This means that rt = ^max 
provides a reasonable estimate for the gravitational tidal 
radius of a galaxy embedded within that subhalo. We 
find that 2% (10%) of our galaxies with L > lOOi© (10 
Lq) have ri/2 > r-max, and that the majority of these 
systems have lost more than 90% of their dark matter 
mass since falling into the host. Physically, these objects 
with ri/2 > ?'max rcprcscnt systems that are losing stellar 
material. We will not explore the observational implica- 
tions of this evaporating population here, but this defini- 
tion may prove useful for future theoretical explorations 
aimed at predicting the fraction of dwarf satellites that 
should be showing signs of ongoing stellar stripping. By 
excluding these stripped galaxies with large values 
from our estimates we are being conservative in the sense 
that we are under-es timating the steal t h popu lation by 
10% for L > IOLq. iPefiarrubia et all (|2008b( ) have ar- 
gued that systems that are losing stellar material such 
as these can actually evolve towards higher mass to light 
ratios once the total bound mass within the stellar core 
decreases by ~ 75%. The central surface brightness typ- 
ically drops by about 2 magnitudes during the process as 
the stellar system expands to adjust for the loss of central 
mass. Further investigation of this possibility and its im- 
plications for the stealth galaxy population is warranted, 
but beyond the scope of this paper. 

Before moving on to our results, we mention that the 
halo finder and the associated definition of halo mass 
used by Mostcr et al. (2009) in our stellar mass assign- 
ment differ slightly from those used in the VL2 catalogs. 
We estimate that this amounts to a ^ 20% difference 
in dark matter halo mass association for any individual 
object, a difference that is not significant given the ex- 
ploratory nature of this work. 

4. RESULTS 



Figure |4] provides a summary of our fiducial model pre- 
dictions compared to current observations: galaxy vs. 
luminosity (left), helio-ccntric distance vs. luminosity 
(middle), and cumulative number vs. luminosity (right). 
The small colored symbols in the left and middle panels 
are model galaxies, with color and symbol type indicating 
three infall redshift bins with Zjnf > 1.5 (red triangles), 
0.5 < Zinf < 1.5 (green diamonds), and z < 0.5 (blue 
plusses). The larger, black squares reproduce the MW 
dSph data from Figure 1. In the middle panel, the helio- 
centric distance for model galaxies is measured from an 
arbitrary point 8 kpc from the host dark matter halo 
center. Our gross results are independent of this choice 
for solar location. 

Model galaxies above the solid line in the left panel of 
Figure m are too diffuse to have been detected. Clearly, 
this population is significant. At fixed luminosity, sys- 
tems above the solid line ~ 30 mag arcsec"^) tend 
to have been accreted more recently (blue plusses were 
accreted since z = 0.5). This trend follows directly from 
our redshift-dependent mapping between L and halo 
mass - at fixed stellar mass, the required halo mass in- 
creases with redshift (Figure 3) and, as discussed above, 
more massive halos tend to host more concentrated stel- 
lar distributions for a given cr* and L. We note that 
for brighter luminosities L > lO^L© our systems tend 
to be smaller in physical size than the observed MW 
dwarfs. As discuss below, this is related to the fact 
that they inhabit more massive dark matter subhalos 
(Afsoo — 3 X 10^) than than appears to be the case for 
the MW dwarfs. The right panel of Figure |4] reveals that 
there is a slight tendency for early-accreted galaxies (red 
triangles with Zinf > 1-5) to be closer to the Sun than 
more recently-accreted galaxies (middle panel) . This bi- 
ases them to host older stellar populations. 

The dotted line in the right panel of Figure S] shows 
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Fig. 5. — Mass within 300 pc as a function of luminosity. Milky 
Way dSph galaxies are shown as black squares. Fiducial model 
galaxies are show as smaller colored points, with the point type 
and color mapped to the time they fell in to the virial radius of the 
VL2 main halo. Earlier accretions are red and recent accretions 
are blue as in Figure 4 and as indicated in the upper panel legend. 
The upper panel presents all predicted satellites and the middle 
panel shows only model satellites that are concentrated enough to 
be detected with current methods, with /i < 30 mag arcsec"^. 
The bottora panel includes only the subset of /.i < 30 mag arcsec"^ 
galaxies that are close enough to have been detected by SDSS, ac- 
cording to the 90% completeness limits in Walsh et al. (2009). The 
M300 masses presented here for MW dSphs are those derived by 
Strigari et al. (2008). We note that subsequent work has suggested 
that t he velocity disper sion (and hence mass) errors on Hercules 
(Herc. Rdeii et a.l.|[2009l) and Wl (B. Willman & M. Geha, private 
communications) are likely underestimated because of membership 
issues. 

the predicted luminosity function of satellites that are 



observable for an SDSS-like survey covering the full sky 
according to the luminosity-distance completeness limits 
of Walsh et al. (2009). This should be compared to the 
data points, which reflect the current M W satellite pop- 
ulatio n corrected for sky-coverage as in iTollerud etaD 
(|2008l ). The uncertainty in the sky-coverage correction 
(associated with the possibility of an anisotropic satellite 
distribution on the sky) is similar to the size of the points 
(ToUerud et al. 2008). We see that the predicted and ob- 
served populations are roughly consistent. The solid line 
shows the predicted luminosity function for all satellites 
within 400 kpc, without any allowance for observational 
incompleteness. The dashed line, on the other hand, 
shows the subset of those galaxies that have peak stellar 
surface densities that are bright enough to be discovered 
with standard techniques. We see that roughly half of 
the systems that are in principle luminous enough to be 
detected with deep surveys like LSST (with L < IOOOLq) 
have peak surface densities that are too diffuse to be seen. 
Specifically, even a survey like LSST, with a very deep 
limiting magnitude, will have difficulties detecting these 
systems without new observing strategies. 

Figure O explores how detection bias affects the Stri- 
gari plot. In the upper panel, we show A/30Q vs. L for 
all model satellites within 400 kpc of the Sun (color/type 
scheme is the same as in Figure S]) compared to the MW 
dSphs (black squares). We see that unlike the data, there 
is a significant population of predicted galaxies with low 
central densities M300 < 5 x IO'^Mq for L < SOOOLq 
. Moreover, while the data follow a nearly common- 
mass relation for M300 vs. L, the model points prefer 
a steeper trend: M300 oc with c ~ 0.15 (as expected 
from abundance matching). The model predictions are 
very similar to those presented in many past CDM-based 
explorations of satellite M300 values (e.g. iBusha et ahl 
I2009t IMufioz et all [20091) . The similari t y betw een our 
model results and those of iMuhoz et all (I2009D . in par- 
ticular, are encouraging. These authors use the same 
VL2 catalog that we use. but they explored a more so- 
phisticated model for assigning light to subhalos. Gener- 
ally, a population oi L < 5000Lq satellite galaxies with 
M3Q0 < 5 X IO^A/q seems to be a fairly robust expec- 
tation for hierarchical models, especially if H2 cooling 
plays a role in the formation star s in surviving galaxy 
halos at 2 = (|Muiioz et al.ll2009[ ). 

The middle panel of Figure [5] includes only those model 
galaxies that have peak surface brightness /i < 30 mag 
arcsec"^. We see that this requirement immediately re- 
moves the population of M300 ^ 5 x IO^Mq objects. The 
lower panel includes only those galaxies that meet both 
the surface brightness requirement and the luminosity- 
distance requirement for SDSS discovery. We see that 
the resultant population of observable model galaxies has 
M^QQ values that arc very much in line with those of the 
known MW dSphs. The model we have adopted therefore 
reproduces both the luminosity function and something 
close to the mass-luminosity trend seen for Milky Way 
dwarfs, once all of the relevant selection bias effects are 
taken into account. 

While we predict that stealth galaxies reside within 
dark matter subhalos that have the smallest total mass 
(as characterized by T4iax or M300), this does not mean 
that these systems will have small dark matter fractions 
within their stellar extents. On the contrary, stealth 



9 



s- 10° 




10^ 10^ 10" 10^ 10^ 10' 
Luminosity [Lq] 




10^ 10^ 10" 10^ 10^ 10' 
Luminosity [Lg] 



Fig. 6. — The dynamical mass-to-light ratio within the 3d half-light radius vs. lu minosity (left) and Vmax vs. luminosity (right). The 
black squares are Milky Way dwarf spheroidal galaxies, with mass estimates are from lWolf et al.] 1 12010) and the Vmax estimates are from 
Wolf et al. (in preparation) . The solid triangles show Fiducial Scenario model galaxies that are concentrated enough to be detected by 
current methods (/i < 30 mag arcsec"'^). Open triangles show stealth galaxies from this same model (^ > 30 mag arcsec^'^). We see that 
stealth galaxies are biased to more dark matter dominated within r-^ii than are their higher surface brightness counterparts. The reason 
is simply that they have larger values. In reality, these model galaxies reside in dark matter halos that are less massive in terms of 
their overall bound mass, Vmax, as can be seen in the right panel. 



galaxies have large half-light radii for their luminosities, 
and therefore sample a large portion of their dark mat- 
ter halos. Indeed stealth galaxies should be more dark 
matter dominated within their stellar radii than higher 
surfaee brightness ultra-faint dwarfs. It is straightfor- 
ward to understand this by considering the integrated 
mass within its 3d half-light radius (Wolf et al. 2010): 
-^(''1/2) = 3 ri/2 crj- For a fixed velocity dispersion 
(which we assume is dictated, on average, by luminosity) 
systems with larger half-light radii will have proportion- 
ally larger half-light masses. Without some theoretical 
prior, this mass by itself tells us little about the total 
dark matter halo mass of the dSph host subhalo, but it 
does enable an easy measure of the mass-to-light ratio 
within the half-light radius T1/2 = ^{^1/2) I {L/^). 

The left panel of Figure H] sh ows T1/2 vs . L for the 
MW dwarfs (open squares from IWolf et all IMol) com- 
pared to our Fiducial model galaxies (filled and open 
triangles). The right panel of Figure [6] shows the same 
simulation data with the same symbol types but now 
with halo Vmax vs. L. The MW dwarf Vmax estimates are 
obtained with a CDM subhalo prior as described in Wolf 
et al. (in preparation). The solid triangles correspond 
to predicted galaxies that are concentrated enough to be 
detected by current methods and open triangles show 
predicted stealth galaxies (/x > 30 mag arcsec" ). As ex- 
pected, the stealth galaxies are biased to more dark mat- 
ter dominated within ri/2 than are their higher surface 
brightness counterparts. The implication is that stealth 
galaxies should have integrated mass-to-light ratios up- 



wards of T 



1/2. 



10 . Again, this is a direct result of the 



fact that they have large half-light radii. As illustrated in 
the right panel, the same galaxies actually reside within 



subhalos that are less massive overall, with Vnax values 
that are systematically smaller than ^ 10 km s^^. It is 
important to remind the reader that the precise values 
of M{rii2) for ri/2 ^ 300 pc in the theory points are 
derived by extrapolating the subhalo profiles inward to 
radii that are not well-resolved in the simulation. On 
the other hand, the Vmax values for the data are extrap- 
olated outward from ri/2- So while neither comparison is 
ideal, the approximations are reasonable, and the qual- 
itative behavior that distinguishes stealth galaxies from 
more readily observable galaxies is expected to be ro- 
bust. Note that the predicted galaxies at the bright end 
of the luminosity distribution are both too small (with 
lower Ti/2 values than the data in the left panel) and too 
massive (higher Vmax values than the data in the right 
panel). This is the same problem discussed above in as- 
sociation with Figure 4, and stems from the fact that the 
VL2 simulation contains several subhalos that are more 
massive than any of the dSph satellites of the Milky Way. 

Of course, if the smallest dark matter halos do not 
contain galaxies at all, then the likelihood for a signif- 
icant stealth galaxy population is much reduced. We 
explore this expectation with our Threshold Scenario, 
which imposes a sharp scale in galaxy formation at 
A^haio = 5 X 1O^M0. Below this scale dark matter ha- 
los are completely devoid of stars (vertical line in Figure 
3). The resultant Strigari plot and luminosity functions 
for this model are shown in Figure [71 Like the Fiducial 
Scenario, the Threshold Scenario also reproduces the ob- 
served satellite luminosity function (right panel). How- 
ever, unlike in the Fiducial Scenario, we now expect only 
a handful of stealth galaxies that remain undiscovered 
(solid vs. dashed lines). The Threshold case also yields 
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Fig. 7. — Threshold Scenario predictions for the Strigari plot (left) and luminosity function (right). The symbols and line types are 
identical to those in Figure 5 and the right panel of Figure 4, respectively. This model imposes a sharp truncation in galaxy formation 
efRciency at A/haio = 5 X 10* A/0, which drastically reduces the expected number of stealth galaxies. 



a Strigari relation that is in reasonable agreement with 
the data (left panel), without appealing to any selec- 
tion bias. Another distinct difference in the Threshold 
Scenario is that all of the low-luminosity galaxies are ex- 
pected to be quite old (or at least to have been accreted 
early Zinf > 1.5). 

One problem with both of our models is that they 
produce too many massive dwarf satellites with A/300 ^ 
2 X 10'' M0 compared to the data (Figures [5l [6l and [7]). 
This has nothing to do with our method of light assign- 
ment. There are simply too many massive subhalos in 
the VL2 halo when compared to the MW satellite popu- 
lation. The same problem is responsible for the fact that 
our model galaxies tend to be too small at i > IO^Lq 
compared to the data in the left panel of Figure 31 These 
model galaxies reside in subhalos that are more massive 
than the subhalos that host the brightest galaxies in the 
Milky Way, and this confines their stars to a character- 
istically smaller radius. Of course, this issue may simply 
reflect cosmic variance in the subhalo populations from 
galaxy to galaxy, but it could be attributed to the fact 
that VL2 represents a slightly more massive halo than 
the Milky Way's dark matter halo. 

It is important to mention that the M300 masses in 
Figures 2, [5] and [7] are taken directly from the spherical 
Jeans modeling of Strigari et al. (2008). A more recent 
analysis of membership in the Hercules dwar f (labeled 
Her in Figures [5] and ^ by lAden et all (120091 ) suggests 
that the actual M300 mass for this system is about a 
factor of ~ 2 lower than the value we have used. Though 
this result does not change the fact that there is no strong 
observed trend between M300 and L in the data, it does 
make the mass of Hercules more difficult to explain in our 
model (and in almost all CDM-base d models to date) . 
In the right panel of Figure 2, the lAden et al] (|2009f ) 
velocity dispersion for Hercules cr^ ~ 3.7kms~^ would 
shift the point at L ~ 3 x lO^L© down to the edge of 
its error-bar ~ clearly low enough that this system would 
not be expected to be stealth according to our definition. 
The fact that the velocity dispersion shifts this much by 
removing or adding a few stars suggests that there is 



a large membership-related systematic error that needs 
to be taken into account in the mass modeling. Q More 
work on the issue of membership and mass modeling in 
this interesting object is certainly warranted. 

5. CONCLUSIONS AND DISCUSSION 

We have argued that there is likely a population of 
low-luminosity satellite galaxies orbiting within the halo 
of the Milky Way that are too diffuse to have been de- 
tected with current star-count surveys, despite the fact 
that they have luminosities similar to those of known 
ultrafaint MW dSphs. These stealth galaxies should 
preferentially inhabit the smallest dark matter subhalos 
that host stars (Knax ^ 15 kms^^). One implication is 
that selection bias (Figures 2 and 5) may play a role ex- 
plaining t he apparent common mass scale for MW dSph 
galaxies (IStrigari et al.l l2008t iPefiarrubia etall l2008al : 
iWalker et al.ll2009HWolf et al.l[2010t ). It also implies that 
searches for the lowest mass 'fossil' galaxies left over 
from reionization may be hindered by surface bright- 
ness limits. This lat ter point was made earlier by 
iBovill fc Ricottil (|2009[ ). According to our fiducial esti- 
mates, potentially half of several hundred satellite galax- 
ies that could be observable by surveys like LSST are 
stealth (Figure 4). 

We developed a plausible estimate for the number and 
character of MW stealth satel lites using the s ubhalo cat- 
alogs of the VL2 simulation (jPiemand et al.l [20081) . We 
assigned light to subhalos by extrapolating the dark mat- 
ter mass-light relations hip required to repr oduce bright 
galaxy number counts (jMoster et al.l [2009^ and we as- 
signed stellar velocity dispersions to each system by 
adopting the empirical relationship between cr^ and L 
for known Milky Way dwarfs. Finally, galaxy sizes were 
computed using the dynamical relationship between 
and CT^ for the measured dark matter halo densities in 
each subhalo (Equation . The resultant model galaxy 
population includes a substantial fraction of ultrafaint 

^ This is not the only systematic uncertainty. As the most elon- 
gated of the dSphs known, Hercules is not particularly well suited 
for spherical Jeans modeling in the first place. 



11 



galaxies that are stealth, with peak surface brightness 
/K > 30 mag arcsec"^. 

We also explored the possibility that there is a sharp 
threshold in galaxy formation at a halo mass of Mhaio = 
5 X lO^MfT-, . This idea fol l ows fr om the common mass con- 
jecture in iStrigari et alJ (j2008| ) and remains viable since 
it reproduces the observed Milky Way satellite luminos- 
ity function as well as a fairly weak M^oo vs. L trend 
without appealing to selection bias (Figure 6). In this 
scenario, all satellite galaxies are born within halos that 
are quite dense, and therefore the number of predicted 
stealth galaxies (which preferentially inhabit the small- 
est dark matter halos) is significantly reduced (Figure 6). 
Moreover, we expect that all of the low-luminosity satel- 
lites will have been accreted since z ~ 1.5 (Figure 6) and 
that they will all host old stellar populations. This is not 
necessarily the case in our Fiducial Scenario, where the 
most distant and low-mass subhalos may host ultrafaint 
galaxies that contain intermediate-age stars. 

A significant shortcoming in our approach is the com- 
plete lack of accounting for the effect of baryons on 
the density structure, distribution, and number counts 
of dwarf galaxy subhalos. The VL2 simulation that 
we have used is purely dissipationless, and therefore 
lacks a central disk component, which would act to en- 
hance subhalo mass loss and select ively eliminate subha- 
los that cross near the disk region (ID'Onghia et al.l[2010l: 
iPenarrubia et al1l201Clf ). At the same time, it also lacks 
the inclusion of any gasdynamical effects such as adia- 
batic contraction, which can incr ease the central densi- 
ties of subhalos with stars. In fact. lWadepuhl fc Springeil 
()2010[) found that subhalo Vmax values tend to increase 
for the most luminous satellites {L > IO^Lq) in cosmo- 
logical simulations. Their simulations self-consistently 
form a central galaxy and produce a satellite population 
that matches the bright end of the luminosity function of 
Milky Way dSphs, so the direction of the effect needs to 
be taken seriously. However, the degree to which baryons 
can enhance the densities of massive subhalos is uncer- 
tain and model specific, as certain wind models can also 
act to evacuate mass f rom the centers of dwarf galaxies 
(|Governato et al.ll2010f ). Ah of these concerns emphasize 
the need for further work in the attempt to make sense 
of observed dwarf galaxy structural properties. At this 
time it is difficult to know how baryons would affect our 
estimates for stealth galaxy counts. For the faint galaxies 
of concern, there is very little baryonic material left to 
have driven contraction, while significant central mass 
loss from winds or global dark matter mass loss from 
tides are both possible. For systems that remain bound, 
these effects will increase the likelihood for them to be 
low surface brightness, as galaxies will tend to expand as 
their global dark matter potential depth decreases. On 
the other hand, mass loss also enhances the rate of dis- 
ruption, and can deplete the overall number of satellites. 
Given that our current model seems to reproduce the 
observed luminosity function (and mass function) fairly 
well, one could argue that depletion has not been too 
much of an issue. 



It is well known that galaxy formation has a primary 
scale - the scale that gives rise to the cutoff in galaxy 
counts at the bright end of the luminosity function. We 
do not know if there is a second scale that operates at 
the low-luminosity end. One implication of our findings 
is that a complete search for very low surface brightness 
satellite galaxies of the Milky Way can help determine 
whether or not there is a second scale in galaxy forma- 
tion. If very few stealth galaxies are discovered, this will 
be an indication that there are no very low mass halos 
that host stars. A similar effect would be seen if there 
were simply a truncation in the power spect rum, as might 
be ex pected in ^ 1 KeV WDM models (IStrigari et al.l 
[2008t iMaccio fc FontanotI [20091 ). In this sense, the dis- 
covery of many stealth galaxies in the halo would provide 
a means to constrain dark matter particle properties in 
a ddition to galaxy f ormat ion physics. 

iKollmeier et al.l (j2Q09( ) have performed kinematic 
follow-up observations of the Pices Overdensity (also 
called Structure J) at a distance of ~ 85 kpc and have 
argued that it may be a very low surface brightness dwarf 
galaxy. If this is true then it represents the first detection 
of a stealth dwarf galaxy in the halo of the Milky Way. 
Unlike the well-known ultrafaint dSphs of the Milky Way, 
which were discovered as overdensities in RGB or MS 
turnoff stars, this system was discovered as an excess 
in RR Lyrae va riables i n the multi-epoch SDSS Stripe- 
82 (jWatkins et al. 2009: iSesar et al.ll2007j ). We caution 
that one potential problem with the stealth-galaxy inter- 
pretation of the Pices Overdensity comes from Sharma 
et al. (in preparation), who have used 2MASS data to 
show that this structure is consistent with being part of 
a larger overdensity of stars, in which case it is unlikely 
to be a bound dwarf. Deeper, wide-field imaging and 
spectroscopic follow-up will be required to determine the 
nature of this interesting structure. 

It is reassuring to note that upcoming deep, time- 
resolved sky surveys and associated follow-up campaigns 
with 30m-class telescopes offer significant hope for the 
discovery of hundreds of new dwarf galaxy companions 
of the Milky Way (ToUerud et al. 2008). Repeated sky 
surveys like Pan-STARRS and LSST may provide the 
best hope for discoveries in the future. Not only will 
they allow the identification of variable tracers, but they 
will also enable concurrent use of bulk proper motions. 
Confident searches within these data offer a means to 
limit the overall count of stealth galaxies that lurk at 
very low surface brightness and to provide unparalleled 
constraints on the efficiency of galaxy formation in the 
smallest dark matter halos. 
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